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ABSTRACT 


Relative  energies  in  kcal/mol 


AEi 

AE2 

E2-E1 

E3-E 

Best  est. 

1.9 

19.8 

14.3 

25.0 

B3LYP 

-8.4 

11.8 

-0.2 

1.9 

M05-2X 

1.4 

16.8 

16.9 

25.4 

It  has  recently  been  pointed  out  that  current  density  functionals  are  inaccurate  for  computing  stereoelectronic  effects 
and  energy  differences  of  isomerization  reactions  and  isodesmic  reactions  involving  alkanes;  this  has  been 
interpreted  as  an  incorrect  prediction  of  medium-range  correlation  energies.  This  letter  shows  that  the  recently 
published  M05-2X  functional  has  good  accuracy  for  all  three  of  the  recently  highlighted  problems,  and  it  should  be 
useful  for  a  wide  variety  of  problems  in  organic  chemistry. 


Density  functional  theory  (DFT)  is  now  the  main  tool  for 
calculating  the  structure  and  energetics  of  complex 
molecular  systems  and  materials.'  Recently,  Grimme 
claimed  that  all  state-of-the-art  density  functionals 
provide  a  qualitatively  incorrect  picture  of  the 
stereoelectronic  effects  in  alkane  isomers,^  and  these 
observations  are  reinforced  by  two  independent  works^’"* 
in  this  journal.  Table  1  compares  Grimme’s  results  for 
four  standard  density  functionals^'^  to  experimentally 
derived^  '"  data  for  the  zero-point  exclusive  energy  of 
«-octane  minus  that  of  2,2,3,3-tetramethylbutane.  All 


four  density  functionals  that  he  studied  get  the  sign 
wrong.  Grimme^  pointed  out  that  a  related  problem  is  the 
increasing  size  of  the  errors  per  bond  in  heats  of 
fomiation  and  alkyl  bond  dissociation  energies  predicted 
by  most  density  functionals  when  molecules  get  larger.^’" 
We  had  pointed  out  earlier  that  the  new  M05-2X  density 
functional  has  much  better  performance  than  all  previous 
density  functionals  for  alkyl  bond  dissociation  energies  as 
a  function  of  alkyl  group  size.'^  Unfortunately,  Grimme 
did  not  test  the  M05-2X  density  functional  for  this 
problem.  In  fact,  we  report  here  that  the  M05-2X 


functional,  with  the  6-311+G(2df,2p)  basis  set,'^  yields 
+2.1  kcal/mol  (+1.4  kcal/mol  with  the  cQZV3P'"'  basis  set 
and  MP2/TZV(d,p)  geometries)  for  the  quantity  in  Table 
1,  which  is  quantitatively  correct. 


of  theory  and  with  the  same  basis  set  as  was  used  for  the  calculation  of 
the  energy.  ^Ref.  12.  -^Present  work  with  the  aug-cc-pVTZ  basis  set.  ^ 
Ref  15. 


Table  1:  Energy  Difference  (kcal/mol)  between  /i-octane  and 


2,2,3,3-tetramethylbutane 

method 

AE  (kcal/mol) 

experiment 

+1.9'^ 

PBE* 

-5.5 

TPSSh 

-6.3'^ 

B3LYP  ® 

-8.4'^ 

BLYP/ 

-9.9^^ 

B3PW91  S 

-7.0  * 

M05-2X ' 

+2.1^ 

M05-2X ' 

+1.4/ 

MP2^ 

+4.6“^ 

^  Refs.  2,  10.  ^  Ref  5.  ^  Calculations  were  performed  with  the 
cQZV3P  basis  set  and  MP2/TZV(d,p)  geometries,  and  results  were 
taken  from  Ref  2.  ^  Ref  6.  ^  Ref  7-9.  ^^Ref  7.  ^  Ref  8.  ^  present 
work  with  the  6-3 1 1+G(2df,2p)  basis  set;  the  geometry  was  optimized  at 
the  same  level  of  theory  and  with  the  same  basis  set  as  was  used  for  the 
calculation  of  the  energy.  *  Ref.  12.  -^Present  work  with  the  cQZV3P 
basis  set  and  MP2/TZV(d,p)  geometries  ^Ref.  15 


Table  2  compares  these  energies  of  reaction  to 
experiment  for  several  common  density  functionals  and 
also  M05-2X.  Clearly,  M05-2X  is  more  accurate  than 
previous  functionals. 

A  related  example,  in  particular  a  case  of  DFT  failing 
to  account  for  stereoelectronic  effects,  was  provided  by 
Schreiner  et  al.,"*  who  compared  the  energies  of  three 
isomers  of  (CH)i2;  see  Figure  1.  Table  3  compares  their 
most  accurate  calculation  and  their  DFT  calculations  to 
our  M05-2X  calculations.  Again,  the  M05-2X  functional 
does  quite  well. 


Wodrich  et  al.^  showed  that  the  systematic  errors  in 
DFT  as  the  alkane  size  is  increased  are  related  to  the 
stabilizing  interaction  of  geminal  methyl  or  methylene 
groups,  an  effect  they  call  “protobranching,”  which  is 
sensitive  to  medium-range  correlation  energy.  A  good 
measure  of  this  effect  is  provided  by  the  energies  of 


reaction  for  reactions  such  as 

n-Ceiiu  +  4  CH4  ^  5  CzHe 

(1) 

or 

n-CgHis  +  6  CH4  ^  7  C2H6 

(2) 

Table  2:  AE  (kcal/mol)  of  reaction  for  isodesmic  reactions  (1) 
and  (2) 


method 

n-hexane 

«-octane 

Experiment 

13.1  “ 

19.8" 

B3LYP  * 

bo 

11.8" 

PBE"^ 

8.9" 

13.9" 

OLYP'^ 

5.9" 

00 

00 

MPWBIK*^ 

9.5 " 

14.4" 

TPSSIKCIS/ 

7.5 " 

11.3  " 

B3PW91^ 

8.0* 

12.0* 

M05-2X  ‘ 

11.5* 

17.2* 

M05-2X ' 

11.2/ 

16.8/ 

MP2^ 

14.1  “ 

Ref  3.  All  DFT  and  MP2  calculation  in  Ref  3  employed  the  aug- 
cc-pVTZ  basis  set.  *  Ref  7-9.  Ref  .5.  ^Ref  16.  ®  Ref  17.  -^Ref  18. 
S  Ref  8.  ^  Present  work  with  the  6-3 1 1+G(2df  2p)  basis  set;  for  all 
calculations  in  this  table  the  geometry  was  optimized  at  the  same  level 


Figure  1  Structures  of  (CH)i2  isomers,  where  1,  2,  and  3 
correpond  to  the  stuctures  of  1,  22,  and  31  in  Ref.  4. 


Table  3.  Energies  (kcal/mol)  of  (CH)^^  isomers  relative  to 
structure  1 


method 

2 

3 

CCSD(T) " 

14.3  * 

25.0* 

BLYP" 

-10.0* 

-11.5* 

G96LYP 

-6.4* 

-6.8* 

KMLYP  " 

41.7* 

B3LYP  / 

-0.2* 

1.9* 

BHandHLYP  S 

7.4* 

14.0* 

B3PW91  * 

14.4* 

19.8* 

B3PW91  * 

15.9' 

22.1  ' 

M05-2x7 

14.0  ‘ 

21.4' 

M05-2x7 

16.9* 

25.4* 

MP2^ 

23.2* 

31.2* 

^  Ref  13,  19,  and  20.  ^  from  Ref  4.  For  each  functional  and  MP2, 
we  give  the  result  with  the  6-3 1 1  +G(d,p)  basis  set.  ^  Ref  7.  ^  Ref  2 1 . 
^Ref  22.‘^Ref  7-9.  ^  Ref  23.  Ref  8.  ^  Present  work  with  the 
6-3 1 1+G(2df,2p)  basis  set;  the  geometry  was  optimized  at  the  same 
level  of  theory  and  with  the  same  basis  set  as  was  used  for  the 
calculation  of  the  energy.  J  Ref  12.  ^  Present  work  with  the 
6-3 1  H-G(d,p)  basis  set;  the  geometry  was  optimized  at  the  same  level  of 
theory  and  with  the  same  basis  set  as  was  used  for  the  calculation  of  the 
energy.  ^Ref.  15. 


We  believe  that  the  sueeess  of  the  M05-2X  funetional 
derives  from  the  design  of  its  funetional  form,*^  building 
primarily  on  work  of  Beeke,^^’^"*  and  from  eonsistent, 
simultaneous  parametrization'^  of  the  M05-2X  exehange 
and  eorrelation  funetionals  against  a  broad  range  of 
aeeurate  data  ineluding  noneovalent  interaetions  for  main 
group  ehemistry.  This  allowed  us  to  better  represent 
medium-range  eorrelation  energy.  We  note  that  the  M05- 
2X  funetional  involves  kinetie  energy  density  in  both  the 
exehange  and  eorrelation  funetionals,  but  unlike  fifth- 
rung  funetionals,^^  it  does  not  involve  terms^^  dependent 
on  virtual  orbitals. 

The  ealeulations  were  earried  out  with  a  modified 
version  of  Gaussian03.^^  M05-2X  is  now  available  in 
release  5.0  of  NWChem.^^  In  addition  it  is  seheduled  to 
be  in  the  next  minor  revision  of  GaussianOS  and  in  the 
next  release  (version  7.0)  of  Jaguar. 

The  reader  is  referred  to  the  original  paper'^  and 
referenees  therein  for  further  details  of  the  M05-2X 
funetional.  In  one  of  our  previous  papers,^^  we  analyzed 
the  noneovalent  interaetion  of  methane  with  benzene  and 
eoneluded  that  the  sueeess  of  M05-2X  for  this 
noneovalent  interaetion  is  due  to  its  improved  eorrelation 
funetional  for  the  deseription  of  medium-range 
eorrelation.  We  noted  in  ref  29  that  although  M05-2X 
does  not  give  the  asymptotie  -Ca/R^  tail  of  the  long-range 
interaetion,  it  agrees  with  CCSD(T)/eomplete -basis 
results  for  CH4-benzene  within  0.17  keal/mol  from  3.4  to 
5.4  A  (the  minimum  is  1.5  keal/mol  at  3.8  A  and  is  4.6 
times  deeper  than  the  value  at  5.4  A). 

We  also  note  that  the  M05-2X  funetional  has  been 
shown  to  outperform  many  other  funetionals  for 
noneovalent  interaetions,  ’  ’  for  torsional  potentials  of 
eonjugated  polyenes,^’  and  for  proton  affinities  of 
eonjugated  polyenes  (a  test  of  its  ability  to  prediet 
polarizabilities  of  eonjugated  double  bonds). It  is  less 
satisfaetory  for  energy  differenees^^  between  eumulenes 
and  poly-ynes,  but  still  better  than  all  other  funetionals  in 
Tables  1-3.^'  M05-2X  also  has  very  high  quantitative 
aeeuraey  for  main-group  thermoehemistry  and  barrier 
heights.'^ 
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